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A digital avioni^.s system referred to as STOLAND has been test-flown 
in the NASA CV-340 aircraft to obtain performance data for time-controlled 
guidance in the manual flight director mode- The advanced system 
components installed in the cockpit included an electronic attitude 
director indicator and an electronic multifunction display. Navigation 
guidance and control computations were all performed on a digital 
computer. A detailed 4D area navigation systems description is given. 

The pilot/system interface and systems operation and performance are also 
described. Approach fllghtpaths were flown which included a 180® turn 
and a 1-min, 5® straight- in approach to 30-m altitude, at which point 
go-around was initiated* Results are presented for 19 approrches. 
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estimated aircraft position In a runway coordinate system 
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A FLIGHT INVESTIGATION OF A 4D AREA NAVIGATION SYSTEM 


CONCEPT FOR STOL AIRCRAFT IN THE TERMINAL AREA 
Frank Neuman, David N. Warner, and Francis J. Moran 
Ames Research Center 


SUMMARY 


A digital avionics system referred to as STOLAND has been test-flown 
in the NASA CV-340 aircraft to obtain performance data for time-controlled 
guidance in the manual flight director mode. The advanced system components 
Installed in thi cockpit included an electronic attitude director indicator 
and an electronic multifunction display. Navigation guidance and control 
computations were all performed on a digital computer. A detailed 4D area 
navigation systems description is given. The pilot/system interface and 
systems operation and performance are also described. Approach fllghtpaths 
were flown which Included a 180® turn and a 1-mln, 5® straight-in approach to 
30-m altitude, at which point go-around was initiated. Results are presented 
for 19 approaches. 


INTRODUCTION 


Studies have shown (e.g., refs, 1 and 2) that short-haul aircraft may 
provide an effective transportation system that can operate into city centers 
and suburban facilities. To provide the detailed data base required for the 
design and development of such a short-haul system, a joint DOT/NASA STOL 
Operating Systems Experiment Program has been initiated. As a part of this 
joint program, Ames Research Center has developed an experiments program with 
the overall objective of providing information that will aid in the choice of 
terminal area guidance, navigation, and control system concepts for short-haul 
aircraft, and in investigating operational procedures. 

In a short-haul transportation system, various levels of avionics systems 
capability may be needed. Simple, low-cost systems may be adequate for 
navigation, guidance, and control of aircraft operating in low-density traffic 
conditions and relatively good weather. More complex and costly automated 
systems may be economically justifiable for operations in hlgh-denslty traffic 
conditions and poor weather. The test data obtained in this program will 
provide a basis for the selection of system capability to meet operational 
requirements (e.g., runway requirements, weather mlnimums, etc.) and will also 
provide means for estimating the system acceptability and system cost. 

A digital avionics system, referred to as STOLAND, has been Installed 
(without servos) in the NASA CV-340 twin-engine transport aircraft. Nineteen 
test flights have been made to obtain preliminary STOLAND performance data 


In the manual flight director mode using time-controlled guidance (l.e., 4D 
area navigation (4D RNAV)). STOLAND Is also Installed (with servos) In the 
powered-llft Augmentor Wing Jet STOL research aircraft (fig. 1) described In 
reference 3 and In a DeHavllland DHC-6 Twin Otter STOL aircraft. Reference 4 
gives an overall functional description of all the STOLAND features for the 
three airplanes In which It was Installed. 

Investigations are being conducted In these aircraft to obtain performance 
data on both simple and sophisticated avionics system concepts. Including auto- 
land, and the corresponding STOL operating procedures. Reference 5 describes 
the technical details of the navigation system and the flight test results 
obtained In the CV-340. The navigation system was specifically designed to 
perform RNAV experiments at Crows Landing NALF, Crows Landing, Calif., which 
has TACAN and an experimental microwave landing system Installed. Refer- 
ence 6 describes the technical details and flight test results of the Ames- 
deslgned 4D RNAV system In the CV-340. This system Included a timed capture 
trajectory and therefore was the first true 4D RNAV system. Reference 7 is a 
brief conference paper on the Sperry-designed 4D RNAV system flight test 
results in the CV-340. Reference 8 ties together the 4D RNAV experience to 
date and also Includes some experience with automatic 4D RNAV in the Augmentor 
Wing Jet STOL research aircraft. 

This report uses the same flight test data as reference 7, but examines 
them In detail. In particular, it distinguishes between left and right 
turning approaches, which brings out the Importance of navigation bias errors 
on fllghtpath distortion. The second contribution of this report is the 
technical description of the 4D RNAV system. 


SYSTEM DESCRIPTION 


STOLAND Is an Integrated digital avionics system which performs all 
terminal area navigation, guidance, and automatic control functions for a 
STOL test vehicle along a curved reference approach fllghtpath. Included In 
the system are the autopilot modes considered standard for commercial trans- 
port aircraft and an autothrottle. The major components of the system are a 
Sperry 1819A general-purpose digital computer and a data adapter that Inter- 
faces all the navigation aids, displays, controls, and servo actuators 
(fig. 2). The navigation aids Include VHF omnirange (VOR), distance measuring 
equipment (DME), tactical air navigation (TACAN), Instrument landing system 
(ILS), microwave modular Instrument landing system described in reference 9 
(HODILS), and a radio altimeter. 

The system components Installed in the cockpit of the aircraft (fig. 3) 
Include a horizontal situation Indicator (HSI), an electronic attitude 
director Indicator (EADI), a multifunction display (MFD), an MFD control 
panel, a mode select panel (MSP), a status panel, and a data entry panel. 
During automatic operation, the pilot monitors the system operation through 
the various cockpit displays. During flight director operation, the pilot 
uses the same set of displays to fly the aircraft along the reference 
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Figure 2.- STOLAND flight-test system. 
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Figure 3-- 3TOLAND cockpit installation. 
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fllghtpath and to monitor the system. A detailed description of the operation i 

of the system is given in reference 4. I 

The navigation system used provides estimates of the aircraft position 
and velocity with respect to a runway coordinate system which has its origin 
at the glide-slope intercept point (fig. 4). The position and velocity 
estimate are generated using ground navigation aid information blended in a 
complementary filter with inertial information obtained from body-mounted I 

accelerometers and attitude sensorn, and air data obtained from a barometric - ] 

altimeter and an airspeed sensor. The ground navigation data are obtained j 

from TACAN except when the aircraft is in MODILS coverage after passing 
point A (fig. 4). The navigation system also estimates wind velocity utiliz- 
ing air data. In the event of a momentary loss of ground radio navigation j 

aid information, navigation is accomplished by dead reckoning using air data. | 

Upon regaining radio information, the system automatically switches back to 
the use of radio data. A detailed description of the navigation system is 
presented in reference 5. 

The RNAV system used during the landing approach is based on a fllghtpath, 
stored in the airborne computer. The fllghtpath is specified by waypoints 
(X,Y,Z coordinates) and associated information such as the radius of turn 
between waypoints and the maximum, minimum, and nominal airspeed between 
waypoints. An Illustration of an approach fllghtpath flown in the CV-340 is 
shown in figure 4. It consists of a long inbound leg (waypoints 1-10), a 180® 
turn to final approach with a 5® glide slope occurring halfway around the turn 
(waypoints 10-12), and a final straight-in approach (waypoints 12-14). A 
brief summary description of the 4D RNAV system is given next. A detailed 
description of the system is given in appendix A. 

The approach guidance is Initiated when the aircraft captures the rear 
extension of the straight line between waypoints (e.g., between waypoints 8 
and 9; see dashed line, fig. 4). At waypoint 8, cont? 'lied time-of-arrlval 
(4D) guidance is initiated. Slightly before waypoint 10, a predictive bank 
angle command is given, and just before waypoint 11, a constant vertical 
acceleration maneuver is performed to acquire the 5® fllghtpath angle. The 
short straight-in section (waypoints 12-13) is the last segment using the 4D 
guidance laws given below. The remaining fllghtpath to flare is flown with 
similar lateral and longitudinal guidance laws except for the system gains, 
which are relatively high from waypoint 13 to flare to assure precise path 
tracking. 

For lateral tracking, the guidance law is: 


K.Y + K,t + (^ 
1 err 2 


( 1 ) 


where (|) equals zero, for a straight line track, and for a circular track, 
P 

V 2 

“ tan~^ — 

P Rg 


(2) 
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For vertical tracking, the guidance law is: 

e “ — h +KY^^^+ — Ih dt (3) 

c V err 4 err v J err 

g g ^ 

As previously stated, 4D RNAV is Initiated at waypoint 8 (fig. A). From 
this point, the system attempts to arrive at waypoint 13 at a given cime. 
Control of arrival time at waypoint 13 is based only on speed control, which 
is provided by controlling the throttle as a function of an airspeed error. 
In the flight director mode, the airspeed command Is displayed on the EADI. 
The airspeed command defined as the algebraic sum of a prescribed 

nominal airspeed (V^j^) and an error that is proportional to an aircraft 
position error (AS): 


V - O.OAAS (m/sec) (A) 

cref AN 

where AS is the distance along the track from the estimated aircraft 
position to a moving target, which represents the desired aircraft position 
(fig. A). As the aircraft arrives at waypoint 8, the target and aircraft 
positions are made to coincide. The computed nominal arrival time at way- 
point 13 is based on the time it would take to fly from waypoint 8 provided 
the aircraft flew the path exactly at the nominal airspeed and provided that 
there was no wind. To account for winds, the position of the moving target 
is recomputed every 10 sec based on the latest estimate of wind velocity 
and direction. This new computed target position assures that the target will 
arrive at waypoint 13 at the nominal arrival time while moving at the nominal 
airspeed. If the wind were changing during the approach, the computed posi- 
tions of the target would have step changes every 10 sec which would result 
in excessive throttle activity. To limit the throttle activity, the time 
rate of change in the value of AS in the above equation is limited to 
6.1 m/sec. For safety and other operational requirements, the airspeed 
command V^^ef limited to preset maximum and minimum airspeeds. 


RESULTS AND DISCUSSION 


As previously noted, the primary purpose of flight tests in the CV-3A0 
was to validate the operation of the STOLAND system and to obtain preliminary 
insight into the performance of ti^e navigation and guidance system. The data 
presented are from a set of 20 simulated IFR (hooded) approaches, which were 
flown at the Crows ^andlng Test Facility (shown in fig. 5 and described in 
derail in ref. 5). It includes a simulated STOL runway, whih is 609.5 m 
long, a CACAN station for navigation, and an experimental microwave landing 
system (MODILS). It also includes a tracking radar to provide tracking data 
for independent postflight calculation of navigation errors. The tracking 
data were smoothed with a minimum mean square second-order curve fit over 
±3 sec to (')btain a best estimate of the actual aircraft position. 


Hi 


Figure 5 


MODILS 



Cround support equipment deployment at Crows Landing. 
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Representative performance of the guidance and navigation systems along 
a typical flightpath will be given. 


Flightpath Description 

Figure 6 shows orthographic projections of (1) the reference approach 
paths and associated waypoints, (2) the aircraft position as measured by 
radar for a typical flight test run, and (3) the position as computed by the 
onboard navigation system. Sixteen left and fou’^ right turning approaches 
were flown on two different days to study the e ects of navigation aid bias 
errors on navigation and guidance. 

The left and right turning flights are presented separately, since the 
errors resulting from the navigation aid bias errors depend on the fl'ghtpath. 
Figure 6(a) shows that the example left approach was initiated to the right 
and above the reference path. Up to and through the turn to the final 
approach, the aircraft remained to the right of the path and then acquired 
the runway centerline, maintaining that course for the remainder of the 
approach. During G/S tracking, the aircraft converged linearly from 30 m to 
20 m above the reference path. This may be attributed to a MODILS elevation 
Mas error of about 0.4®. The major error prior to MODILS acquisition can be 
attributed to the effect of the TACAN DME bias. Figure 6(b) shows similar 
results for a right turning approach. 


Wind Environment 

Following is a description of the x and y wind components which were 
encoun ered during the test flights, as estimated by the navigation system. 
For the first day, the x wind along the runway was a relatively steady tail 
wind of 7.3 m/sec. The y wind was -6 m/sec (from the left) at the 506-m 
approach altitude and reduced almost linearly to 0 m/sec at the altitude 
at which the approach was terminated (30 m). On the second day, the x wind 
component was ^ :ail wind of 10 m/sec, and the y wind component linearly 
changed with decreasing altitude from -4.6 m/sec to 1.5 m/sec at the approach 
termination altitude. 


Navigation and Guidance Errors 

The navigation errors (shown in figs. 7. 9, and 11) include the ground 
navaid and airborne receiver signal errors, off nominal atmosphere effects on 
the altimeter, errors in the radar tracking data, and the basic navigation 
system errors resulting from software/hardware mechanization. The correspond- 
ing guidance errors are showri in figures 8, 10, and 12. In these figures, the 
abscissa is the distance (S) along the track. For orientation, the waypoints 
are labeled on the figures. Waypoints 10, 11, and 12 describe the final turn, 
and waypoint 13 is the final waypoint for 4D RNAV. Data are given for a typi- 
cal flight along with the envelope of all data. The derivation of the equa- 
tions for the navigation error calculations is given in appendix B. Navigation 
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errors are positive if the estimated position along the flightpath is ahead 
or to the left of the actual (radar-derived) position. Guidance errors are 
positive if the estimated position la ahead or to the left of the reference 
position. 

Cvoastvaak navigation evvova- The crosstrack navigation error is the 
difference between the crosstrack positions as measured by the precision radar 
and the STOLAND navigation estimate. For the left-turning flightpaths of 
figure 6(a), the envelope of lateral navigation errors is shown in figure 7(a). 
At the initiation of the approach at waypoint 8, errors as large as 200 m 
occur. These errors converge to less than 70 m at the initiation of the turn 
at waypoint 10, where they start to Increase again to values as large as 150 m. 
Examination of the data indicate that these navigation errors result from 
TACAN bias errors in both range and a2imuth« A short time after passing 
waypoint 10, the transition from TACAN to MODILS navigation is initiated. 
Navigation errors then converge smoothly to less than 15 m after transition 
to MODILS is completed. 

For the right-turning flightpaths of figure 6(b), the navigation errors 
show a different, but equally consistent, pattern (fig. 7(b)). Here, the 
TACAN DME error causes the aircraft to be flown on the left side of the 
nominal reference path. The errors after transition from TACAN to MODILS 
are eventually reduced to values similar to those shown in figure 7(a), since 
the flightpaths are common from waypoint 12 on. It will be noted that the 
envelope of figure 7(b) is narrower than that of figure 7(a). This is simply 
because the 4 right-turning flights were conducted on one day, whereas 16 left- 
turning flights were conducted on two different days; experience has shown that 
the navalds have different bias errors on different days. In contrast, a 
comparison of the data for the four left and four right approaches on the same 
day showed that they both had narrow envelopes. 

Cvoaatvaok guidance errors^ The error shown in figure 8(a) is the 
difference between the onboard estimate o^ ositlon and the minimum distance 
to the reference flightpath. The envelope ux the cross track guidance error 
shows large errors at the initiation of the approach due to the use of 
different capture maneuvers. As a result of guidance activity, these errors 
converge to smaller values before switching to MODILS. To explain what 
happens between waypoints 10 and 11, one must consider the effect of transi- 
tioning between navalds with different bias errors. After switching to the 
more accurate navigation aid shortly after waypoint 10, the new position 
estimate rather quickly Indicates that the aircraft is not on the reference 
path but to the right of it. This accounts for the sudden Increase in guidance 
error, for which the pilot has to correct. The low-gain 4D RNAV lateral 
guidance law slowly guides the aircraft back to the reference path until, 
after passing waypoint 11, the guidance errors converge to small values. 
Finally, the high-gain localizer tracking law, which takes over after way- 
point 13, converges the guidance errors to a narrow range of ±20 m when the 
aircraft is 1600 m from touchdown. Figure 8(b) shows a trend similar to fig- 
ure 8(a). Positive guidance errors in this case mean that the onboard 
navigation system estimates that the aircraft slightly on the left of the 
reference path. 
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Vertioal navigation errors- The envelope of the vertical navigation 
errors, figure 9(a), shows errors as large as 24 m at initiation of the 
approach* The vertical navigation errors are always positive and are probably 
the result of a bias in the baro-altimeter. It should be noted that the baro*- 
altlmeter reference was set before each approach, based on information radioed 
from the control wower. This setting only assures that the barometric alti- 
tude coincides with the geometric altitude at the runway level* As a result, 
different bias errors occurred on each approach, which makes the envelope in 
figure 9(a) fairly wide* After transition to MODILS and the start of the 
descent at waypoint 11, the baro-altlmeter measurement is slowly combined with, 
and, after 1 min, fully replaced by the MODILS data, to prevent a step change 
in estimated altitude at the initiation of glide-slope tracking. Due to a 
constant 0.4® MODILS elevation bias error, ^ the vertical navigation error 
past waypoint 13 converges linearly with distance from the MODILS elevation 
antenna. After waypoint 13, the width of the error envelope is approximately 
8 m* Figure 9(b) illustrates the same sequence of events for right turns: 

(1) there was an initial barometric bias error of about 10 m past waypoint 10, 

(2) during the 1-min transition to MODILS data, the errors increased due to 
the MODILS elevation bias error, and (3) fur the straight-ln tracking past 
waypoint 13, the errors decreased linearly as before* A comparison of the 
data for the four left and four right approaches on the same day showed that 
the barometric altitudes had the same bias errors. 

Vertical guidance errors- The envelope of vertical guidance errors 
(fig. lO(.a)^ shows a range of errors as large as 15 m at the initiation of 
the approach at waypoint 8 and is generally above the desired path* The 
magnitude of the error represented by the envelope width remains approxi- 
mately constant between waypoints 8 and 10. The vertical guidance error in 
figure 10 is blanked out during the -5.0® flightpath angle capture maneuver, 
which starts approximately at waypoint 11. The vertical guidance error during 
this maneuver has no meaning since the capture law between path segments does 
not use the altitude error, but it is a closed loop law based on a linearly 
varying flightpath angle command, 0^. = K(y - Between waypoints 13 and 

14, the vertical guidance error envelope converges to about ±3 m as a result 
of the high-gain guidance law and high-gain navigation filters used during the 
final straight-ln approach. The guidance errors in figure 10(b) are very 
similar to those in figure 10(a), although the character is somewhat different 
since the corresponding navigation errors have different spatial changes. 

Alongtraak navigation errors- The alongtrack navigation errors, shown 
in figure 11, contribute to time-of-arrival errors at the final waypoint. 

As expected, the navigation errors have different characteristics for left 
and right approaches, since the bias errors distort the measured flightpaths 
differently. The error trend is even different in the early common portion 
of the flightpath (past waypoint 13), although the MODILS azimuth bias errors 


^In later flights, the MODILS elevation bias error was reduced consider- 
ably. 

^In an advanced version of this system, the capture law was replaced by 
guidance along a circular arc as in equation (3) where h^^^ is the devia- 
tion from the arc. 
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Figure 11.- Alongtrack navigation errors. 
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did not change. The difference Is explained by the long time constant of 
the X-navlgatlon filter, and the fact that the filter had different bias 
error Inputs before reaching waypoint 13. 

Alongtraak guidance errova~ Typical alongtrack guidance errors are shown 
In figure 12 for left and right turns. These were used to command airspeed 
corrections on the EAOI based on the control law V(. ■ .V^ + 0.04 ASl. Again, 
the characteristics of the curves are different for left and right turns, 
and similar for turns In the same direction. However, the dependence of the 
guidance errors on the navigation errors and the pilot control was not analyzed 
In detail, because (1) some of the wind calculations for turns, which were 
used for nominal tlme-of-arrlval calculations, were later found to be in error 
(see eqs. (A41) and (A42) In appendix A); (2) only the rate-limited AS,(ASj^), 
has been recorded, and (3) the throttle position, which Is the pilot response 
to the Indicated error, was not recorded. 

Comparison of figures 12(a) and 12(b) shows that the average time error 
at the final waypoint (expressed as AS^ errors) Is smaller for the right 
turns than for the left turns. A reason for this can be seen from figures 6(a) 
and 6(b), which are typical of all approaches. The left turning fllghtpath 
(fig. 6(a)) Is several hundred meters longer than the reference fllghtpath. 

This difference must be made up by speed control, but the final waypoint (13) 

Is too close for a complete correction due to the limited speed variations 
available. In contrast, for the right-turning fllghtpaths, the aircraft flew 
a shorter path on the Inside of the reference path, and, therefore, It did not 
have to make up for Increased path length with speed control. 


Tiroe-of -Arrival Errors at Waypoint 13 

Figure 13 Is a histogram of the tlme-of -arrival errors at waypoint 13 for 
the simulated Instrument (hooded) approaches. For these tests, the mean 
time-of-arrlval error is 3.1 sec (late) with 2o deviation of ±4.3 sec. The 
mean tlme-of-arrlval error obtained ..urlng these tests may result from the 
TACAN range error which caused the actual longitudinal distance flown to be 
longer titan the reference path. Additional data are required to establish 
the system performance for all TACAN and MODILS error combinations. 

It Is Interesting to note that current manual guidance techniques enable 
air traffic controllers to deliver CTOL aircraft to the runway within about 
±15 sec of the predicted arrival time (ref. 10). This capability corresponds 
to p single runway acceptance rate of about 40 IFR arrivals per hour using 
current separation standards. Using the Improved capability of the automatic 
tlme-of-arrlval guidance system described here. It would be possible to 
Increase the runway acceptance rate by about 40 percent (see ref. 10). 
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CONCLUSIONS 


Results are presented for 20 flight director approaches made during an 
Investigation of a 4D RNAV system using the NASA CV-340 aircraft. Results 
of these limited tests led to the following conclusions: 

1. Blended radio /Inertial navigation using TACAN and a microwave 
scanning beam landing guidance system (MODUS) permitted a smooth transition 
from area navigation (TACAN) to precision terminal navigation (MODILS). 

2. Time of arrival at a point about 2 n. ml. (3.70 km) from touch- 
down was about 3 sec ±4 sec (2o) later than the computed nominal arrival time. 

3. Navigation system bias errors (rather than system noise) are 
primarily responsible for the distortion of the fllghtpaths and the cha. r 
of the guidance errors. 



APPENDIX A 


STOUND 4D RNAV 


INTRODUCTION 

The purpose of this appendix Is to give a technical overview description 
of the STOLAND 4D RNAV system which was used to obtain the experimental data 
In this report. The experimental tests hsa limited objectives, and the 
system was designed with these objectives in mind. The primary purpose was 
Co test path guidance concepts, Including :lme-of-arrlval control, once the 
reference path was acquired. The system v/ss designed primarily for automatic 
guidance. However, for the tests described In the body of the report, flight 
director commands were obtained through further processing of the guidance 
system commands. For completeness, this processing will also be described In 
the appendix. 

The navigation system, which computes estimates of position x, y, z, 

ground speeds x, y, z, and wind components }L, y^ In a runway coordinate 
system, has been described elsewhere (ref. 6). In this appc.ndlx, equctt^ns 
will be given that transform the navigation quantities into values relative 
to the flightpath, as needed In the guidance computations. 


FLIGHTPATH DESIGN CONSIDERATIONS 


The reference flightpath must be capable of being flown under all 
environmental conditions, particularly wind and turbulence, given the 
aircraft performance limitations and control authority limits. In order to 
limit the complexity of In-f light computations, a few requirements are made 
on the geometry of the trajectory and on the location of the waypoints. 

Segments between waypoints must be either completely linear or 
completely circular. Circular flightpath segments must have large enough 
radii so that the roll angle command will not exceed some specified limit, 
say 30”, for wind conditions likely to be encountered with the maximum 
specified airspeed (for tlme-of-arrival control), allowing a margin for 
path error correction. 

Stralght-llne and circular flightpath sections must be connected in such 
a manner that the circles are tangent to the straight lines. The first 
circular segment must be entered from a straight line segment. After the 
first segment, there is no restriction on the number of subsequent .circular 
segments that may be contiguous except that there must be no discontinuities 
in direction of the trajectory between a circular segment and any contiguous 
segment. Segments between waypoints must represent flightpaths that the 
aircraft can fly within the specified airspeeds. Stralght-llne flightpath 
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segments may follov each ocher without a connecting circular segment. In this 
rase, the system will perform appropriate capture maneuvers between line seg- 
ments. In Che present STOIAND configuration, the final line segment of the 
4D RNAV flightpath must be an extension of the glide slope and localizer tra > 
to be flown to touchdown. The prescribed airspeed and airspeed limits vary 
linearly between the reference values defined at the waypoints. 


Flightpath Definition 

Pctpametera that define the flightpath- In STOLAND, provision is made for 
four prestored fllghtpaths, each defined by a maximum of 30 waypoints. Each 
waypoint is specified by seven quantities (x, y, z, R, ^min* ^ref^* 

The X, y, and z are coordinates of Che waypoint with respect to a terminal 
area (runway orlentc - ) coordinate frame. The radius of curvature R Is 
specified when a specific flightpath curvature Is desired. The polarity of R 
determines the direction of the turn>. positive for right turn and negative 
for left turn. The waypoint at which R is specified represents the start 
of that specific curved path. If zero curvature Is desired (R ^ »), then 
the notation used to recognize that a straight-line segment is desired Is to 
enter zero in the R column. The terms V^gx* '^mln* '^ref specify the 
maximum and minimum speeds allowable and the reference speed at the waypoints. 

Auxiliary parametera- The Information entered In the computer defines 
the flightpath, but it Is insufficient to actually fly it and to display the 
flightpath on the multifunction display. For these purposes. Information 
from two sequential waypoints Is used to calculate auxiliary quantities, 
namely, the glide-slope angles, heading, and direction sine and cosine for 
straight-line sections, and tae coordinates of the center of turn, heading 
change anglr in the turn, and other quantities for circular line sections. 

The auxiliary parameters used for the 3D and 4D RNAV are computed in 
an initialization routine following waypoint entry, after one of the four 
fllghtpaths is selected for display. If waypoints are changed in flight or 
on the ground, a reinitialization is performed. 

The following parameters are computed for straight-line segments 
(fig. Al) . The distance d^ is the distance along the specified flightpath 
between waypoints m and n: 


d 

n 


^ n m n 



(Al) 


The next three parameters arc constants that are related to the orientation 
of a line segment. The course is 



(A2) 
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Direction sine and cosine are 


■ y*. 

n ID 


- sin ij;. 


b - - Q 08 Cj 

n n 


(A3) 

(A4) 


The following parameters are computed for curved sef;ments (fig. A2). 
Coordinates of the center of a circular segment are obtained from the 

radius of segment n, Ym, the coordinates at the start of segment n, 
and the course at point (x^, Yjn)- If is entered positive for right 

turn and negative for left turn, the equations are 


■'b ■ “n 
n 

'’O. ■ "n 


(A5) 


(A6) 


For circular flightpaths, the course at the end of the turn the angle 

of turn Ai|/, and the arc lengt!i between waypoints are defined below. 


tb tan 
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(A7) 


For a right turn (R positive), 
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For a left turn (R negative). 


= - ip -2tt if ip - 4f >0 
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(A8a) 


(A8b) 


d = R At|i 
n ' n ' 


(A9) 


The nominal flightpath angle for both turning and straight-line segments is 

(-V^) 
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Another auxiliary parameter is the velocity gradient. Ii is used 

to calculate the nominal speeds and speed boundaries between waypoints. 


D 

n 


V - V 
ref ref 
n m 

d 

n 


(All) 


Note that the acceleration, dv/dt, is not constant in this case. Similar 
expressions hold between the minimum and maximum speed boundaries. 


3D Fllghtpath Guidance Computations 

The ooject of the fllghtpath guidance is to fly closely to the reference 
path, in spite of navigation errors and physical disturbances, without 
excessive control activity. Figure A3 shows the final portion of a typical 
fllghtpath used in the flight tests. It defines the guidance error quanti- 
ties. It ends with a 180° turn onto the final leg. The last 90® of this final 
turn descends at a glide-slope angle of 5®. The final waypoint is about 1 min 
from touchdown. 


Independent of path geometry, the path tracking guidance laws remain 
essentially the same, except for the reference quantities of fllghtpath angle 
and nominal roll angle (f»p and some gains. 

For lateral tracking, the guidance law uses the crosstrack error Yerr 
to null the error and the crosstrack velocity Ye^r lateral fllghtpath 

damping: 


= Kj + K V 
rerr 2-^err 


(A12) 


where the predicted roll angle is 

((» = 0 
P 

for a straight-line track, and 



(A13) 


(A14) 


for a circular track, where the ground speed is calculated from navigation 
data 


V 

g 


= + ^2 


(A15) 


and g *=■ 10 m/aac^. 
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For vertical control, the pitch guidance lav uses the altitude error 
*err null the error and the Inertial fllghtpath angle error Yerr 
vertical fllghtpath damping Is 


e 

c 


V *err 
8 


+ KuY 
n'err 



where the fllghtpath angle error is 


(A16) 


’err ' ’ret ‘ ’l 
n 


(A17) 


The inertial flightpath angle is calculated from navigation data. For 

small angles » 

A 

yi~f- (A 18 ) 

g 


The remaining error quantities in the above guidance equations, Yerr> 
Yerr’ ^err> computed from fllghtpath and navigation data as follows. 
The aircraft path errors Yerr ^err computed as horizontal and 

vertical distances from the fllghtpath as defined by a plane perpendicular 
to the fllghtpath which contains the center of gravity of the airplane. For 
a straight-'llne reference fllghtpath, the lateral deviation is (fig. A4) 


y 


err 


a (x - X ) - b (y - y ) 
n n n n 


(A19) 


which is positive when the aircraft is on the left of the reference fllghtpath. 

The crosstrack velocity is calculated from the navigation system ground 
speed estimates x and y. For a straight-line path 


err 


a X - b y 
n n 


(A20) 


From figure A5, the lateral deviation from a circular fllghtpath is 

Ygrr * slgn(R) ~ ^)^ + ~ y)^ - |R|] (A21) 

defined so that, in agreement with the straight-line fllghtpath, the error 
is positive if the aircraft is on the left of the fllghtpath for either 
right or left turns. The crosstrack speed error for a circular path is 


err 


(x_ - i)x + (y 
slgn(R) „ = 

yj(x - i)2 + (y 


R 
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y)y 

W 


(A22) 
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The vertical deviation from the reference fllghtpath is (fig. A6) 
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err 
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(A23) 


Both straight and circular fllghtpaths require drng« the range to the next 
waypoint » for this error calculation. For a straight-line fllghtpath segment, 
the range to the next waypoint Is 


or 
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> b 
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(A24) 


(A25) 


provided that the aircraft Is close to the fllghtpath. That these are 
reasonable approximations for small crosstrack errors can be seen from the 
exact expressions which include the cross track distance (fig* A7) 
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(A27) 


For a circular fllghtpath, the range from the aircraft to the next waypoint 
is (fig. A8) 


d « (A28) 
rng ' ' 

The complete calculations of equation (A28) are equivalent to those of^ 
equations (A7) to (A9) except that Xj^ and yj^ are replaced by x and y, 
and is replaced by 


Path Segment Capture 

Since the aircraft responds at a finite rate, and since passenger 
comfort must be observed, transition maneuvers are required between fllghtpath 
segments. Guidance law switching is accomplished by monitoring the distance 
drng* Before entering the next path segment, a capture maneuver is performed 
which will result in a smooth entry to that segment. This is important when 
the fllghtpath angle or turn radius Is changed, and when transitioning between 
straight-line segments of different courses* 
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Figure A7.- Range to next waypoint geometry for 
straight-line segment. 


32 





Hovizontal capture- In order to lead a turn between segments, the bank 
angle connnand Is Issued somewhat before the waypoint is reached. The decision 
to switch between horizontal path segments Is based upon the best estimate of 
ground speed and range to the waypoint according to the relationship 


KV > d , K - 2.9 sec (A29) 

g rng 

where K was found experimentally to result in a smooth transition. 

Vertical capture- The system switched between vertical path elements 
when the predicted pitch maneuver will result In a nominal normal acceleration 
during the transition. From figure A9 the switching range Is 


and further. 


rng 




(A30) 


(A31) 


where R Is the circular pitch maneuver radius, the acceleration 
nacc “ 0.69 m/sec^, Vg is the ground speed, and Ay Is the difference 
between fllghtpath angles, which Is also equal to the predicted change in 
pitch angle during the maneuver. The transition In vertical guidance and the 
predicted pitch maneuver are switched when 

n > 0.69 m/sec^ 
acc - 


or 


V 2ay 
g - 

2 X 0.69 


> d 

rng 


(A32) 


From the above discussion, it is clear that vertical and circular 
captures are each Initiated at different distances from the next 

waypoint. Therefore, fllghtpath design for the present system is limited 
to having either a change of fllghtpath angle or a change of fllghtpath radius 
at a waypoint. To complete the vertical capture maneuver, the distance d^j^g 
for changing glide-slope angle must be shorter than the distances between the 
last and the present wr.vpoint and between the present and the next waypoint 
(fig. A8) . The distance d^j^g for roll Initiation is usually relatively small 
compared to that for glide-slope angle changes. To remove the above restric- 
tions would result in a much more complex logic structure of the program. This 
concludes the description of the 3D path tracking guidance. Tlme-of-arrlval 
control, the fourth dimension of the guidance problem will now be discussed. 
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(a) Right turn. 


(b) Left turn. 


Figure A8.- Range to next waypoint geometry 
for circular flightpath. 



# Present FPA 
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Figure A9.~ Vertical capture maneuver. 
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Tlme-of-Arrlval Control 

The object of tlme-of-arVlval control is to permit the aircraft to arrive 
at a specified waypoint within a closely defined time interval, in order to 
permit hlgh*density sequencing operations for aircraft in the terminal area. 

The tlme-of-arrival control at the final waypoint is based on speed control 

only. Path stretching capability has not been provided. Speed control for 

4D RNAV is achieved by providing an airspeed command to the aircraft speed 

control system. | 

I 

The airspeed command is defined to be the sum of a prescribed I 

nominal airspeed and an error which is proportional to an aircraft | 

position error | 

I 

V i: “ V^„ + K AS , K - 0.04 sec'^ (A33) j 

cref AN 6 ’ 6 I 

The term AS is the distance along the track from the aircraft position to 
a moving phantom target position (fig. A3), which is positive when the phantom 
is ahead of the aircraft. The reference airspeed V^^ef calculated via 
equation (A33) is also governed to fall between specified speed limits which 
are dictated by aircraft structural design safety margin and also by ATC 
cons ider a t ions « 

The generated phantom target position is the desired aircraft position. 

Initially, at the aircraft’s arrival at the first selected waypoint, the 
phantom and aircraft positions coincide. The computed time it will take from 
the selected capture waypoint to the final waypoint is based on the aircraft 
flying the complete trajectory exactly at the nominal airspeed, provided the 
wind does not change. The arrival time is thus fixed when the aircraft passes 
the capture waypoint, and the ^Ime is displayed to the pilot on the MFD. For 
test purposes a time delay or advance can be added to the nominal flight time 
before engaging the 4D system by entering the number of seconds delay as 

WTE » on the STOLAND keyboard and then reinitializing the system by 

entering WDE » 1. This simulates late or early arrival at the first selected 
waypoint, which the system will try to make up by speed control. 

From the foregoing discussion it is apparent that the nominal time of 
arrival will not be known until passing the first waypoint, nor can it be ‘ 

preassigned. In an operational system this would not ''ul^ill the requirements. 

One solution, investigated at Ames Research Center, employs a timed, continu- 
ously recomputed, capture flightpath yCef- 9) to a selected waypoint so as to 
achieve the desired time of arrival. 

Every 10 sec the position of the phantom aircraft is recomputed based on 
the latest estimate of wind velocity and direction, such that the critical 
waypoint arrival time will be satisfied with the aircraft flying with the 
leasf. deviation from the nominal airspeed Hence, the system attempts 

to meet the critical waypoint arrival time flying at the nominal airspeed. J. 

As long as the aircraft and the desired target position coincide, the aircraft. \ 

will deviate from nominal airspeed only due to subsequent changes in the " 

estimated wind velocity. | 

I 

i 

A" 
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Calculation of the Distance Between Phantom and Aircraft 


The oomponente of LS- Figure 10 Illustrates the notation for the 
distance components employed In computing the position error AS 


AS » d. + S . + S „ 


(A34) 


The aircraft position Is measured with respect to the waypoint next ahead 
of the aircraft, denoted as 1^+1, and the phantom target position is 
measured with respect to the waypoint xj which Is irmedlately behind the 
target. Distances are measured in the direction of flight along the track. 
Desired arrival times {t^} at the Intermediate waypoints are recomputed every 
10 sec, based on the nominal airspeed the latest estimate of wind 

velocity and distances between waypoints {d^}. The generated phantom 

target position will therefore make discrete Jumps resulting from these wind 
changes. However, it represents the best estimate of where the aircraft 
should be to satisfy the critical arrival time, assuming the wind to remain 
constant thereafter. To prevent such fast Jumps, which would cause excessive 
throttle activity, AS is rate limited to 6.1 m/sec. In addition, limits 
on position and rate of change of the airspeed command signal ^cref 
to reduce Jumps in the throttle command signal. 

Computation of phantom target position within a path segmey\t^ The 
specification of the phantom target position will now be outlined. The time 
at which the phantom enters a fllghtpath segment is determined by the target 
time computations given in a following section. 

The problem of generating the phantom trajectory is essentially that of 
venerating a nominal ground speed Vq^ along a specified fllghtpath for a 
gxven wind estimate. The nominal airspeed is specified for each waypoint 
so that the ground speed profile Is the vector summation of the wind with the 
airspeed profile, where the airsjeed profile between waypoints Includes any 
nominal acceleration programs. We are interested in the solution of the 
following 'Integral equation: 



The integrand is ground speed, which is a function of position along the tra- 
jectory, which in turn is a function of time t*. The nominal ground speed 

Is equal to the vector sum of the nominal airspeed and the wind estimate: 


GN 




- sln^ a + V cos a 


(A36) 


The geometrical relationship of equation (A36) Is lllustraced in figure 11. 
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For the computation of the phantom trajectory. It la assumed that 


so that 


V 2 sln^ a << V. 2 
w AN 


= V.„ + V cos a 
GN AN w 


(A37) 


(A38) 


(which assumes a relatively small crab angle). 

The nominal airspeed Van can be expressed as a linear function of 
the distance along the segment i 


• ''aS.1 * 


(A39) 


where is the prescribed nominal air velocity at waypoint i, and 

is a velocity gradient that defines a nominal acceleration program, 
equation (All). Similar expressions apply to the maximum and minimum speed 
limits as illustrated in figure A12. To accommodate a changing wind angle a 
along a curved segment, let 

Si 

a(Si) = cii - ^ (A40) 


where is the angle of with respect to the direction of the flight- 

path segment immediately beyond waypoint (i), is the distance along the 
segment, and is the radius of curvature of the segment with the proper 

sign for right or left turn. Substituting equation (A40) into (A38) and 
expanding gives 


S S 

1 1 


(A41) 


which is valid for both straight and curved segments (fig. A13) - The terms 
and are the crosstrack and the alongtrack wind components at waypoint 
i, respectively: 


M s V sin a. = -a.V + b.V 
i w i I wx 1 wy 


(A42) 


N. “ cos 01 . = b.V + a.V 
i w i i wx 1 wy 


(A43) 


where v^^ and Vyy are the filtered x and y runway axes components of 
V^,, the wind estimate. 


^In the CV-340 implementation! the sign of R was not considered. 
^In the CV-340, a^ instead of -3;^ was used. 
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The terms and are the direction sine b cosine of the waypoint- 
specified trajectory (eqs* (A3) and (A4)). The expression for VQj^(Si) 
given in equation (A41) could be used as the integrand of equation (A5), and 
the integration performed continuously in real time would yield the desired 
position function However, if we are willing to accept a phantom 

trajectory which satisfies arrival time and endpoint velocities but which may 
not follow the precise profile given in equation (A35), a simpler function 
for employed which has the advantage of avoiding excessive 

computational complexity. 


=4.T<^' " (''ak.1^ + - \) * - \Y ^ - hj <*“> 


where 


y AN,i^ i^ AN,i^“^l i^+l/ 


It At, 2 
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(A45) 


and 


“2d j V AN , im Irp ANji^i"! i_'^l/ 


T T 




At4 


(A46) 


The constants A^.^ and are calculated from incorporating the endpoint 

constraints on the time and velocity in equation (A44). The parameter 
At^, the time between waypoints i and i + 1 on the nominal groundspeed 
profile, is determined in the following section. 

For the special case where the segment is linear and is constant 

over the segment, Aj^ = = 0, and the target position becomes simply 



Since some of the parameters used in the above equations (Nj^, t^, Atj^) are 
recomputed at the 10-sec wind-estimate update cycle, Sj^ will generally 

make discrete jumps at these times. Methods to reduce the undesirable 
effects on throttle activity arising from this computational technique are 
discussed in reference 10. 


Arrival Time Computation 

The final waypoint is designated as ’’critical** in that the aircraft is 
required to arrive there at the time which was calculated on initial path 
acquisition. Desired ’’arrival times’* for the waypoints leading to the 
critical waypoint in order to make good this requirement are computed every 
10 sec. 
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Arrival times {tj^} for the waypoints leading to the critical waypoint 
J are obtained in terms of computed segment times and the critical 

arrival time tj as follows: 

J-1 

t » t - T At.,i < J (A48) 

i ‘J jti J 

Arrival times are required only for the waypoints in the range 


1^ < i < J 


(A49) 


See figure A14 for an illustration of the computation. 

The time required to fly from waypoint i to waypoint i + 1 at a 
nominal ground velocity Vqj^ along the trajectory is 

_d. 


ds. 

"'i J 


(A50) 


where is the distance variable measured along the segment from waypoint 

i, and d^ is the distance between waypoints i and 1+1. The term Vqj^ was 
defined in the last section as the vector sum of the prescribed nominal air 
velocity and the estimated wind velocity (eq. (A38)): 


V 


GN 


+ V cos a 
AN w 


where a is the angle of V with respect to Vqj| as Illustrated in fig- 
ure All. ^ 

A further Justification for using the approximate expression (A38) for 
^GN equation (A50) is that small errors in the computation of Atj^ do 
not affect the time of arrival at the critical waypoint as long as it is 
possible to meet the arrival time within the airspeed limits. Small errors 
in computing Atj^ result in small deviations in ground speed from the 
nominal value. This has been shown to result in acceptable tolerances on 
time-of~arrlval control. 

For cases other than linear segments with constant the solution of 

equation (A50) is approximated by rectangular egration. The expression for 
^GN^^i) ^8 derived in the last section is: 


S S 

■ ''an.i ^ Vl * "l r 7 ^ "l r7 


(A51) 
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For case shown 14 through tg must be 
computed t^O = J = critical time 



9 

t4=tio- E ^ti 
j = 4 


Figure A14.- Illustration of computation of t^. 
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where ^ Is a data input nominal airspeed, Is a velocity gradient 

(eq. (All)), and and are computed at 10-sec Intervals from updated 
wind estimates and constants obtained In the Initialization routine 
(eqs. (A42) and (A43)). Therefore 


'■Jo 


ds. 


cos (S^/R^) + Mj sin (S^/R^) 


(A52) 


For the special case where the segment Is linear and Is constant. 


Ati = v 


AN,i ^ '‘l 


(A53) 


and this simplification Is always employed in this case. Otherwise, rectangu- 
lar Integration will be employed to approximate equation (AS2) as follows: 


At, 


AS, 


K-l 
k=i I 




(dj - KAS^) - 




(A54) 


where Is a judiciously chosen integration increment, and K is the 

greatest integer less than dj/ASj^. For a circular segment, AS^ is 
selected to be the distance traveled in a 15® turn: 


AS. = ~ (A55) 

By thus constraining kASj^/Rj^ to take on fixed discrete values, the sine and 
cosine functions in VQj^(kAS^) can be stored in very brief look-up tables. 

For a linear segment (with changing V^), ASj^ is selected to be 152.4 m. 

If the aircraft is greatly behind schedule, it is either accelerating 
or flying at maximum permissible airspeed attempting to catch the target 
position. Similarly, if it is greatly ahead of schedule, it is either 
decelerating or flying at minimum permitted airspeed. The amount by which 
the aircraft is expected to be in error from the nominal arrival time due to 
the imposed maximum and minimum speed limits is calculated as follows for 
presentation on the multifunction display. 

The minimum time required to get to critical waypoint J is given by 


min 

A,J 


min 


+ 


J-1 



A 



(A56) 
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where 
tin 


^mln 


is the minimum time required to go to the next waypoint, and 


Atf*“ is the minimum time required to go from waypoint 


i to waypoint i + 1. 


An analogous expression can be written for the maximum time possible in going 
to waypoint J, denoted by tJ®j. The computation of At^^" and At^** is 

identical to the computation of At^ described above, except that Vy^(S^) 

is replaced with Vjf«*(S£) and respectively. Both velocity 

functions are required to vary linearly between waypoints and thus can be 
expressed in terms of their values at the waypoints as for V^(Sj). 


by 


The minimum time required to get to the next waypoint, 1^+1, is given 


min 




ds. 


A,i^+l 




(A57) 


where i® the distance to waypoint ly^+1 at the 10-sec computation 

Interval, and V^**(S^) la the same as VGjj(Si) except that V^(Sj^) Is 
replaced with V^*(Si). If segment 1^ Is linear and 


vr*(Si) 


min 




max 

''a 1 i 
’A A 


Is constant. 


(A58) 


In general » 


min 

^A.l^+l 
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‘A 


ds. 






^ax ^ + N. cos(S,/R, ) + M sln(S/R. ) 

All* 1* 1 1a 1 ^A ^ 


(A59) 


which Is also approximated by rectangular Integration as in equation (A54), 

and t?®* . Is obtained by an analogous expression. 

A. 1^+1 

The aircraft is expected to be able to arrive at critical waypoint J 
at the prescribed arrival time tj If the remaining time tj - t lies 
within the range 


min 

*^A,J 


< (t 


t) < t 


max 

A,J 


(A60) 
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where t is the current real time. Therefore, if 

then the aircraft is estimated to arrive late by 

*^late ” ~ *^A,J 


Similarly, if 


^ ^ ^max 

*^A,J 


Then the aircraft is estimated to arrive early by 


t = t - t'"^^ - t 

early J A,J 


(A61) 


(A62) 


(A63) 


(A64) 


The Flight Director System 

The 4D RNAV system is designed for automatic operation with the pitch, 
roll, and speed commands driving control servos. The control servos were not 
installed, however, on the first test aircraft which provided the data in 
this report. For this reason, tiie 4D RNAV commands were converted to flight 
director commands. To complete the system description for the data reported 
here, the conditioning of the control law signals, described earlier and 
suitable for input to the flight director system, will be described next. 

The pitch command, computed by the automatic control law (eq. (A16)) 
is first limited as an inverse function of airspeed (see fig. A15) • The 
limited pitch command is then compared with the actual pitch of the airplane, 
which is washed out over 16 sec to a] low the aircraft to assume new trim 
conditions when fllghtpath angle or airspeed is changed. The pitch rate 
sensor signal is added to the vertical reference pitch signal through a gain 
in order to provide damping in tracking the pitch command bar. To prevent the 
flight director bar from jumping when the system is engaged, the pitch angle 
is synchronized whenever a new guidance mode is engaged. 

As shown in figure Alb, the roll command output from the automatic 
control system is limited to 30®, filtered with a 1-sec time constant low- 
pass filter, and limited to a 6® per second roll rate to further reduce the 
effects of the navigation noise. The filtered roll command is then compared 
with the measured roll angle, and roll rate information is added for damping. 
The resulting roll error is limited and scaled for suitable presentation on 
the display. 
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For tlme*of-arrival control without autothrottle, a speed control flight 
director Is required. Therefore, an error signal In the form of a black bar 
was displayed on the EADI, which represents the rate-limited difference 
between commanded and measured airspeed. For the conventional aircraft, the 
CV-340, this was sufficient for reasonable tlme-of-arrlval control. 
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APPENDIX B 


NAVIGATION ERRORS IN MOVING RECTANGULAR COORDINATES 


To study the interplay of navigation and guidance errors, it is convenient 
to express the navigation errors in the Lame coordinate system as guiuance 
errors* Since guidance errors are calculated as alongtrack and crosstrack 
errors, the same convention is used for expressing navigation errors in this 
report: 

The following definitions are needed (see fig. Bl): 

positive if estimated position ahead of radar 

positive if estimated position to right of radar 

positive if estimated altitude above radar 

unit vectors defining radar coordinate frame 

unit vectors defining moving coordinate frame 

components of reference flightpath position measured 
in radar coordinate frame 

components of estimated airplane position measured in 
radar coordinate frame 

Derivation of the navigation error is as follows: 


^NE 

"ne 

^R* ^R* "r 

X, y, z 
*R* ^R’ ^R 
X, y, 2 



(Bl) 


N.E. 


= N.E. I 

Vr"r 


(B2) 


Equation (B2) states that the vector navigation error in the radar 
frame of reference is equal to the vector navigation error in the rotating 
frame. Expansion of equation (B2) gives 


(X - + (y - y^)y^ + 


(B3) 
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( 0 ,- 8204 ) 


(- 11626 ,- 8204 ) 


Figure Bl.- Navigation error calculation geometry 



Substitution of equation (Bl) Into (B3) gives 

^NE^ ” (x - Xj^)(si(iy + cij;x) + (y - yj^) (cij^y - stpx) + (z . 


Therefore, 




Xj^)cij; - (y 
Xj^)sip ■>• (y 






- Z 


z)z 

(B4) 


(B5) 


Equation (B5) requires the angle ip. For r.he initial straight part of the 
flightpath, it is 180®, For the final straight portion of the path, it is 
0®. For the intermediate semicircular portion, it may be found as follows: 


(«„ - -f 

Vcxj - + (y^ - y^)2 


(B6) 


Substitution of y from equation (Bl) into equation (B6) gives 


ctpy + siix. 


^ 

^ V(Xr - Xc>2 + (y^ - y^)2 


(B7) 


R -^c^ 


Therefore, 


sip = 


\ \ 


y/^\ - + (yR - y,)2 


(B8a) 


or 


ci/; = 




V(Xr - X ^)2 + (y^ - y ^)2 


tan Ip 



(B8b) 


(B9) 


Equation (B8) gives values for and cij^ directly hut requires a 
square root, whereas equation (B9) requires an ATAN and then taking the sine 
and cosine. The results aie tabulated on the following page. 
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